TABLE 1 

General Characteristics of Sample Bursts. 



Burst 


Trigger 


Detector 


Model 


No. 


Time 


Average Fluence 


Name 3, 


No. 


No. b 


Used c 


Spectra 


Interval (s) 


-Speak (keV) (erg cm _2 ) d 



3B 910421 
3B 910425 
3B 910503 
3B 910522 
3B 910601 
3B 910619 
3B 910627 
3B 910717 
3B 910807 
3B 910814C 
3B 910814 
3B 911031 
3B 911118 

4(3 




3ken PoweP LcJvv^ Suiuk E 



8.2e-6 
4.8e-5 
7.9e-5 
3.5e-5 
2.1e-5 
4.1e-5 
1.5e-5 
8.3e-6 
2.6e-5 
3.0e-5 
7.8e-5 
3.0e-5 
5.6e-5 
2.0e-5 
2.4e-5 
4.2e-5 
1.8e-5 
4.1e-5 
1.3e-5 
8.3e-5 
5.7e-6 
l.le-4 
5.8e-5 
6.4c-5 
1.0e-4 
2.4e-5 
l.le-4 
9.9e-6 
2.6e-5 
4.5e-5 
1.3e-5 
6.3e-5 
4.4e-5 
2.0e-5 
2.2e-5 
5.7e-5 
4.9e-5 
1.2e-5 
2.5e-5 
1.9e-5 
3.7e-5 
1.2e-4 
2.4e-5 
2.4e-5 
1.2e-4 
7.1e-5 
1.6e-5 



Log. Probabilities 
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TABLE 1— Continued 



Burst 


Trigger 


Detector 
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No. 
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Average 
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a Burst names arc from Meegan ct al. 1996, and continue with a '4B' prefix for bursts that occurred after 
the end of the 3B Catalog. 

b Entries with multiple detectors indicate that summed data were used for the analysis. Entries prefixed 
by 'S' indicate that SD data were used instead of missing LAD data. 

C 'GRB': the spectral form of Band ct al. (1993); 'BPL': broken power law; and 'COMP': power law with 
a high— energy exponential cut— off, in the same -E pe ak parametrization as the GRB model. 

d Fluences as determined over the selected time interval for the given spectral form; differs from the 3B 
Catalog values. 
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TABLE 2 

Analysis of High Energy Spectral Index. 
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a Burst names are the same as in Tabic 1. Bursts fitted with the COMP model are excluded 
from this table. 

b Probability for \ 2 , resulting from a fit to each of the two models for the behavior of (3 with 

time. Very small or zero values indicate that the model was unacceptable. 
c Given as a signed probability. Please sec the text for a full description. 
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ABSTRACT 

We report on the temporal behavior of the high-energy power law continuum 
component of gamma-ray burst spectra with data obtained by the Burst and 
Transient Source Experiment. We have selected 126 high fluence and high flux 
bursts from the beginning of the mission up until the present. Much of the data 
were obtained with the Large Area Detectors, which have nearly all-sky coverage, 
excellent sensitivity over two decades of energy and moderate energy resolution, 
ideal for continuum spectra studies of a large sample of bursts at high time 
resolution. At least 8 spectra from each burst were fitted with a spectral form 
that consisted of a low-energy power law, a spectral break at middle energies and 
a high-energy continuum. In most bursts (122), the high-energy continuum was 
consistent with a power law. The evolution of the fitted high-energy power-law 
index over the selected spectra for each burst is inconsistent with a constant for 
34% of the total sample. The sample distribution of the average value for the 
index from each burst is fairly narrow, centered on —2.12. A linear trend in time 
is ruled out for only 20% of the bursts, with hard-to-soft evolution dominating 
the sample (100 events). The distribution for the total change in the power-law 
index over the duration of a burst peaks at the value —0.37, and is characterized 
by a median absolute deviation of 0.39, arguing that a single physical process 
is involved. We present analyses of the correlation of the power-law index with 
time, burst intensity and low-energy time evolution. In general, we confirm the 
general hard-to-soft spectral evolution observed in the low-energy component 
of the continuum, while presenting evidence that this evolution is different in 
nature from that of the rest of the continuum. 

Subject headings: gamma-rays: bursts 
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1. Introduction 

In the first six years of operation, the Burst and Transient Source Experiment 
(BATSE), on board the Compton Gamma Ray Observatory (CGRO), has accumulated 
a vast amount of spectral data on gamma-ray bursts. Although the BATSE Large Area 
Detectors (LADs) have only moderate energy resolution compared with the Spectroscopy 
Detectors (SDs), they have unprecedented effective area over their entire energy range (28 
keV -1.8 MeV). By studying spectroscopy data from the LADs for bright events such as 
those reported on by [Ford et al. (1995)] , who used SD data, we can track the evolution of 



fitted spectral parameters with finer time resolution, and we can extend the analysis to 
fainter events. In this paper, we analyze 126 bursts at high time resolution, with more than 
8 spectra per event, concentrating on the higher-energy behavior, where it was difficult for 
the SDs to obtain good statistics. 

As with much of the field of GRB studies, theoretical modeling of continuum spectral 
emission naturally breaks into two periods: before and after the publication of the first 
BATSE results flMeegan et al. 1992|) . The paired observation of burst isotropy on the 
sky along with an inhomogeneous distribution of events with brightness, and presumably 
distance, has established the conclusion that GRBs occur much farther away, and are 
consequently much brighter, than previously expected. Instead of comprising a nearby 
Galactic disk population, burst sources either reside in a very large Galactic halo or else 
they are truly cosmological (we will not consider here another possible scenario: that bursts 
may arise in a local heliospheric halo, such as the Oort cloud [e.g.: |Bickert fc Greiner 1994 ; 



but also see: |Glarke, Blaes, fc Tremaine 1994|| ). Such an uncertainty in distance has had 



dire theoretical consequences; no single model has surfaced that can accommodate both 
distance scales, since such a model would have to account for luminousities that differ by 
~ 10 orders of magnitude. The early theoretical work was dominated by the physics of 
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strong-magnetic field Galactic-disk neutron stars (see |Harding 1991 , for a review), which 



has as its basis the efficient mechanism of quantum synchrotron emission. Of course, 
energization of these systems was a crucial problem, in that the emission timescales are on 
the order of 10~ 17 s, for a typical field strength of 10 12 G required to produce a cyclotron 
absorption line fundamental at ~ 20 keV, as observed in X-ray pulsars (|Voges et al. 1982] ). 
Nevertheless, continuum modeling of then-current spectral data enjoyed a moderate success 
(many references in [Ho, Epstein fc Fenimore 1992j ) . 



All this began to break down with the placement of burst sources no closer than a large 
Galactic halo, as most of the strong-field models have restrictive luminosity constraints. 
Cosmological burst emission scenarios proposed to date are less predictive, but have 
had little time yet to mature. For the most part, interest has been focused on merging 
neutron stars, since the total energy budget is about right for very distant events. What 
happens after the merger is what distinguishes the models from each other. A simple 
fireball was proposed by many workers flCavallo fc Rees 1978] , Goodman 1986| , |Piran 19"94|) . 



Non-thermal emission, such as is observed in GRB spectra, is very difficult to produce 
in an optically-thick source, although as a fireball expands and becomes optically thin, a 
high-energy power-law component becomes possible. However, it was soon realized that 
in the environment of two colliding compact objects, baryon contamination of the fireball 
would pose a problem, diverting energy from the direct production of fireball radiation 
into the acceleration of material (see discussion in [Fishman fc Meegan 1995|) . In order 



to address this problem, several workers proposed that the observed gamma-ray emission 
originates not in the original event but is a by-product of the kinetic energy gained at 
the expense of the fireball. Maximal acceleration of the explosion products leads to a 
relativistic blast front, which can cause shocks when colliding with interstellar material, 
either by encountering dense knots or eventually by sweeping up matter in the path of the 
shock front ( [Rees fc Meszaros 1992] , [Meszaros fc Rees 1993| ). Shocks can also arise internal 
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to the outgoing relativistic wind, in the case where the central engine is variable ([Rees fc 



Meszaros 1994j , [Paczyriski fc Xu 1994| ). It is important to note that the energy distribution 



of the shock-accelerated particles that gives rise to the observed emission is not predicted in 
any of these models; however, the distributions can be inferred from observation. The most 
efficient radiation mechanism is synchrotron, which produces a characteristic low-energy 
power-law behavior flKatz 19"9~""| , [1'avani 1996| ). The high-energy spectral shape for this 



model comes from the distribution of Lorentz factors for the baryons arising in the shock, 
which is typically a broken power law. Dispersion of blast-front velocities will give rise to 
observable hard-to-soft spectral evolution, both in individual pulses, as well as over the 
course of the entire event. Some of this behavior has been noted by [Ford et al. (1995) 



however, the opposite behavior is also seen, as well as a mixture of both. 

Apart from the details of individual theoretical models, what new can be learned from 
analysis of spectral data? First, we have the well-known observation that GRB spectra are 
non-thermal. There is good evidence that some time-averaged GRB spectra are composed 
of power-law emission to several 10s of MeV in energy ( [Matz et al. 1985| , [Hanlon et al.| 



1995). Burst emission indeed reaches very high energies, as evidenced by the single 18 GeV 



photon observed by EGRET ( [Hurley et al. 1 9*9*3 ), albeit at a considerable delay from the 



initial outburst. This alone can say much about the distribution of particles doing the 
emitting, as well as the possible optical depth. Other than a multi-temperature blackbody, 
which can mimic a power-law spectrum over a limited energy range, non-thermal emission 
arises from non-thermal particles. The evolution of the particle distribution, by cooling, for 
example, bears a simple relationship to the evolution of the emission for many radiation 
mechanisms. In the fireball model, optical depths are much greater than unity during the 
phase in which the matter gets accelerated. Thermal emission from the fireball is not 
observed in the gamma-ray band (although it may be visible in X-rays below ~ 20 keV, 
Preece et al. 19961) . In any cosmological model, it is very difficult to avoid conditions that 
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will rapidly lead to large optical depths via the photon-photon pair production process. 
This occurs in the collision of two photons where the product of their energies is greater 
than 2ml/(l — cost/)), ip being the angle between the photons' directions and m e = 511 keV 
is the rest mass of the electron. Many bursts have substantial emission at 500 keV and 
greater, so if the high-energy emission is not to be quenched by a runaway pair-fireball, the 
emission must be highly beamed. The high energy power-law index and its time evolution 
should constrain the mechanism through which particles are giving up their energy in 
emission, as well as reflecting the behavior of the injection mechanism. Cases where the 
high-energy component comes and goes within a burst or is absent altogether may represent 
quenching by a mechanism that rapidly increases the optical depth, such as photon-photon 
pair production. In this case, it is expected that the intensity should drop during periods 
of quenched emission, or in other words, there should be a hardness-intensity correlation. 

In this paper, we will present a study of the time-evolution of burst spectra, 
concentrating on the high-energy power-law component. In §2 we discuss the burst sample 
selection and the details of the spectral fitting analysis. The results are covered in §3 
and their implications are discussed in §4. In the Appendices, we summarize the general 
characteristics of BATSE and then discuss in detail the energy calibration procedure for the 
LADs, which has made the current work possible. 

2. Analysis Methodology 

In order to have a sample of bursts with at least 8 spectra with count rates high 
enough to obtain well-determined spectral parameters, we selected a subset of bright bursts 
based upon either the total fluence or peak flux, as determined from the LAD 4-channel 
discriminator data flMeegan et al. 1996| ). We required a fluence (> 20 keV) greater than 
4 x 10 -5 erg cm -2 . However, the set of bursts for which the fluence can be calculated 
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is limited by several considerations, such as data availability, telemetry gaps in the data 
coverage and possible contamination of portions of some bursts with other active sources (in 
particular, with solar flares in the first year of the mission). Thus, we made an additional 
selection of those bursts which had a peak flux from 50 - 300 keV on the 256 ms timescale 
in the 3B Catalog (and later) above 10 photon s _1 cm~ 2 . Each burst was then binned in 
time, so that each spectrum to be analyzed had a signal to noise ratio (SNR) of at least 45 
in the typically 28 to 1800 keV energy range of the HERB data (High Energy Resolution 
Burst data: for a description of the instrument and spectroscopy datatypes, please see 
Appendix A). Bursts with less than 8 spectra after binning were dropped from the sample. 
Most spectra in bright bursts are well in excess of this SNR, which guarantees > 2a of 
signal per energy resolution element, assuming a flat count spectrum. Roughly 20 resolution 
elements (= 5E, the FWHM of the detector energy resolution) are required to cover the 
typical LAD energy range, thus the 128-channel HERB spectra are over-resolved in energy. 
LAD data types other than HERB are under-resolved, which is why HERB is preferred for 
spectroscopy. Some bursts did not have complete coverage in the HERB data (especially 
before a flight software revision that allowed longer accumulations during quiescent portions 
of a burst), in which case we used other available data, as discussed below. There are 126 
BATSE bursts in our sample matching these criteria. 

Background was determined independently for each channel, typically using spectra 
from within ±1000 s of the burst trigger (giving at least three background HER spectra 
before and after the burst). The form of the background model was a fourth-order 
polynomial in each energy channel, where the fitted rates are time-averaged over each 
spectral accumulation, rather than determined at the centers. This was done to avoid 
underestimating the background rate at a peak or overestimating it at a valley. The SNR 
was determined by comparison with the chosen background model, interpolated to the time 
of the accumulated spectrum. 
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Spectra were fitted by one of several spectral forms, depending upon the best fit 
obtained to the average spectrum over the entire burst. The primary spectral form we used 
is the function of [Band et al. 1992| (GRB, in table |l|), which consists of two smoothly-joined 



power-laws: 

/(S) = A(S/100) a exp(-S(2 + a)/S peaj£ ) 

if E < (a-(3)E peak /(2 + a), (1) 

and f(E) = A{(a-(3)E peak /[100(2 + a)]}^eiq>((3-a)(E/100f 

if E > (a-{3)E peak /(2 + a). 

The two power-law indices, a and (3, are constrained such that the resulting model is always 
concave downwards (a > j3; our definition includes a possible minus sign for each index). 
If, in addition, the high-energy power-law index (/?) is less than —2, the model peaks in 
vT v (that is, E 2 times the photon spectrum) within the BATSE energy range. The model 
is parameterized by the energy of the peak in vT v (Speak), rather than the energy of the 
break between the power laws (S = E pe&k (2 + a) /{a — /?)). If the fitted value of (3 is 
very negative, roughly less than —5, the spectral form approaches that of unsaturated 
inverse- Compton thermal emission (|Rybicki fc Lightman 1979|) , a low-energy power-law 
with an exponential cut-off (COMP, in table |l|). This can be viewed as a generalization of 
the spectral form of optically-thin thermal bremsstrahlung (neglecting any Gaunt factor), 
which has a low-energy power law index of —1. The GRB spectral form is a continuous 
function that does not allow a sharp spectral break, so that in cases where Speak (< So) is 
close to the high end of the energy range for the data, (3 may not be well-determined. For 
such cases, we used instead a simple broken power law (BPL, in table [I|), in order to force 
Speak = S , usually resulting in acceptable fits to the high-energy component. 

Since we are concerned in this paper about the high-energy power-law behavior, we 
made a number of tests to be sure that our choices of spectral models do not affect the end 
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result. To do this, we fit several trial bursts with several different models and compared 
the resulting fits. The simplest test of the robustness of our fitting procedure was to fit a 
single power law to each spectrum above a fixed cut-off energy that was determined by 
the maximum over the entire set of fitted spectra in the burst of the value of the break 
energy E between the two power-law components. This eliminated any affect the fit to 
the low-energy data might have on the fitted value of /3. That is, curvature in the global 
model fit may tend to pull the local fit of the high-energy power-law index to a larger or 
smaller value, depending on how well the actual data tolerate the curvature. For example, 
the data may break more sharply than the model, which leads to a fitted value for f3 that is 
steeper than it should be. Conversely, in the broken power-law model, with no curvature 
built in, the high energy power law index may be pulled to a shallower value than the data 
require. Generally, the resulting time-history of the fitted parameters are consistent to 
within one-sigma errors. However, some differences were apparent when we compared the 
time-histories of the fitted high-energy power-law component between these two models, 
when both were applied to the same burst, as can be seen in figure [I]. The average values 
of the fitted power-law indices (weighted by the errors) over the entire burst were slightly 
different (/9 av e = —2.25, for the GRB model fit; = —2.16, for the BPL), while the underlying 
pattern of the time-history of the parameters were similar. So while the time evolution of 
the high energy portion of the spectrum could be reliably traced by the fitted parameter 
for each model, there remains some ambiguity in the average high-energy slope. This effect 
should be worse for larger average values of -E pea k: the curvature inherent in the GRB 
model tends to restrict the range of energies available for determining (3. The broken 
power-law model is plagued by a different problem: with an energy resolution (FWHM) 
of approximately 20% at 511 keV, we usually cannot determine the exact position of the 
break energy using LAD count spectral data. 

With the fitted values of the break energy and (3 possibly closely correlated, the 
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reported la error on each parameter is only part of the story. That is, the errors are 
most accurately determined from a multi-dimensional x 2 contour plot for the correlated 
parameters, as seen in figure || The contours represent Ax 2 values appropriate for one 
parameter of interest, so that the la contour is at Ax 2 = 1 (for this figure only; usually, 
one would be interested in both parameters jointly, resulting in larger contours). The ID 
la error limits are formed by the maximum and minimum of the error ellipse projected 
onto the axis of the parameter being considered. The actual la errors reported here are 
obtained from the diagonal elements of the covariance matrix for each fit; this is equivalent 
to Ax 2 = 1, with the additional assumption that the fitted parameter value lies in the 
center of the error interval. By taking into account the joint error between the parameters, 
Ax 2 is increased to 2.3, so that the fitted values of the high-energy power law index can be 
reconciled to within one or two sigma between the two different spectral forms. 

The fact that we obtain acceptable fits with different spectral models reflects on the 
ambiguity of the forward-folding process. Given the detector response, a count recorded 
in a given data bin could have come from a photon of any number of different energies, 
all greater than or approximately equal to the nominal energy range of the data bin. 
The dominant component of the response at low energies is the resolution-broadened 
photo-peak, centered on the photon energy. On top of this are counts derived from 
incomplete absorption of higher-energy photons in the detector, the off-diagonal component 
of the response. Consistent with the constraints imposed by the detector model, including 
especially the energy resolution, a given photon model folded through the detector response 
matrix will redistribute the predicted counts to best agree with the observed data. Thus, 
the solution to the forward-folding spectral fitting problem is not unique. 

Table [I] summarizes global aspects of the fits performed for each burst. We use the 
3B catalog name ( pVlccgan et al. 1996 ) and BATSE trigger number to identify each burst, 
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followed by the number of the detector with the smallest zenith angle with respect to the 
source, the spectral model used for fitting, the number of fitted spectra, the time interval 
selected for fitting, the average of the fitted values for -E pca k and the fluence, summed over 
the fitted spectra. In cases where there are two or more detectors reported in the third 
column, a summed 16 energy channel data type (MER) was used, usually for lengthy 
events which ran out of HERB memory before the end of the burst. For a small number of 
cases where other data types were absent, we use SD 256 energy channel data (SHERB); 
these are indicated in column three with an 'S' appended before the detector number. The 
three models used in our analyses are indicated by their respective mnemonics (introduced 
above) in column four. The COMP spectral form has one less parameter than the others: 
there is no fitted high-energy power-law index. However, each of the models shares three 
corresponding parameters: amplitude, low-energy power-law index and -E pca k (or spectral 
break energy for the broken power-law model). In the last two columns we indicate the 
average value for Speak in keV and the total fluence for the fitted interval in erg cm~ 2 . 
Notice that three of the four bursts that required the COMP model did so because the high 
energy power law was completely unconstrained; indeed, for these bursts -E pca k was also 
unconstrained, as the average value is far greater than the energy of the highest channel 
available in the data (typically 1800 keV). In the following analyses, we shall exclude these 
four bursts, since no trend in the high-energy power law index can be determined with our 
data. 

3. Observations 

We should like to know several things concerning the behavior of the high-energy 
power-law as a function of time. First of all, is it constant? If not, does the index change 
smoothly with time, as with the hard-to-soft spectral evolution observed in the E pea k 
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parameter by |Ford et al. (1995)|? If the behavior is not smooth, is it correlated with other 



observable features in the burst time history, such as the instantaneous flux or the evolution 
of the low-energy spectral parameters? To investigate these questions, we subjected the 
fitted values of the high-energy power-law index to several statistical tests, and evaluated 
the probability the outcome of each could have occurred randomly. The results of our 
analyses, shown in table 0, are described below. Each row of the table is indexed in the first 
column by the trigger name from table |I|. For each burst, this is followed by the weighted 
average of /3, the probability that a constant (3 describes the data, the probability that a 
linear trend in (3 describes the data, the slope from a linear fit to the time series of (3, and 
the probabilities that the fitted values of (3 are correlated with time, the burst time history 
or with the time series of -Ep ea k- 

To start with, we would like to test the hypothesis that (3 is a constant over the entire 
burst. In order to do this, we first computed a weighted average of the fitted values of the 
high-energy power-law index (which we will denote as (3, regardless of which model we used 
for the fit) over the time interval selected for each burst. The weight applied to each term 
in the average is the square of the 1-sigma error, o"j, of the fit: 



(Ji J I i \cr. 



In cases where the fit resulted in an undetermined value for (3 for an individual spectrum, 
the value was thrown out of the weighted average. It should be noted that, with weighting 
of the individual values, as well as the elimination of undetermined values, the result is 
different from the value of (3 obtained from a fit to the integrated spectrum. The third 
column in table ^] gives the probability for \ 2 obtained by subtracting the weighted average 
from the actual fitted values in each burst. The x 2_va l ues are calculated assuming the 
model, and thus a small value (such as < 10 -4 ) indicates a problem with the assumption 
and thus the likelihood that the model is false. A histogram of the logarithm of these 
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probabilities in figure |3] (dotted line), shows that for some bursts, at least, a constant (3 is 
consistent. What is not shown are the 30 bursts for which the probability is essentially 
zero. Including the bursts for which the log. of the probability is less than —4, we have 
42 out of 122, or 34% of the total sample, that are not consistent with a simple, constant 
model in (3. It is extremely unlikely that this distribution occurs randomly. 

The distribution of (3 ave , shown as a histogram in figure £|, improves on earlier work 
by |Band et al. 1993] , with a larger sample and better statistics per burst. However, the 
resulting values from these two studies cannot be compared directly, since here we have 
weighted each fitted value of (3 by the parameter error, while in the previous study the 
fits were made to average spectra, which are implicitly weighted by intensity. Finally, the 
sample sets are different: the selection of events in [Band et al. 1993| was based upon peak 
counts, not fluence or peak flux, since these were unknown at the time. The median value 
for the sample is —2.12, with an absolute deviation width of wadcv = j§ X/fLi \ Xj — x \ = 
0.23 (where x me d represents the median, which minimizes the absolute deviation), compared 
with the standard deviation of 0.30. The distribution has an extended negative tail that 
gives it a skew value of —0.73 (the skew is defined as the dimensionless third moment of the 
distribution, and is for a Gaussian), large compared with the expected standard deviation 
of the skew of yl5/iV = 0.35 for a purely Gaussian distribution. Given the large variation 
of other spectral parameters, such as E pca ^ which has a distribution at least as wide as 
the range of possible values, it is surprising that the high-energy behavior is so restricted. 
Plotted over the total distribution in figure ^ is a histogram of those bursts for which (3 is 
consistent with being constant (log. probability > —4 from figure |3|). 

Obviously, a constant value of (3 is not acceptable for many bursts. A clear example 
of this is presented in figure [|, which shows the time history of (3 during 3B911118 and 
is an example of general hard-to-soft spectral evolution in (3. The Spearman rank-order 
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correlation of f3 with time is given in column 6 of table ^. The correlation coefficient r is 
distributed between —1 and 1, and can be converted through the combination 

to a Student's t-distribution for N — 2 degrees of freedom. Unlike the x 2 probabilities, 
correlation coefficients that are not consistent with roughly a normal distribution around 
reject the null hypothesis that no correlation exists; therefore, small probabilities indicate 
significant correlation. The probabilities associated with r, calculated using equation 3 
along with the number of spectra fitted (JV) from column 5 of table [I], reveal that a trend 
in the data exists for at least 21 of the events at the 10~ 3 significance level or smaller. This 
is a robust estimator for correlation; it indicates when a correlation is almost certainly 
present. However, the Spearman test does not take into account the errors for each point, 
so if there are a large number of outliers with large errors in the sample, the test will come 
up with poor results. Figure presents the distribution of the time correlation coefficients 
(solid line). The bulk of the distribution consists of negative correlations, indicating an 
anti-correlation of the power-law index with time, or hard-to-soft spectral evolution. 

A linear fit to the time history of /3 also indicates whether there is a monotonic trend 
in the data, while accurately treating the errors in the fitted power-law indices. The fifth 
column of table |2| gives the linear coefficient, or slope, of such a fit, having the units of 
change in f3 per unit time, or s _1 . The sign is such that hard-to-soft spectral evolution (f3 
grows more negative in time) results in a negative slope. The x 2 probability for this fit 
is given in the fourth column of the table and the distribution is also plotted on figure [| 
(solid line). In 24 cases out of the total sample, the log. probability was less than —4, 
indicating that the linear trend was a poor model of the data for those events. Comparing 
this result to that for the model of constant f3, however, more bursts had acceptable fits to 
a linear trend at the same significance level (98 compared with 80 out of 122). There are far 
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more cases of hard-to-soft spectral evolution (100) than there are for soft-to-hard evolution, 
which was already evident in figure |6|. The first spectrum in many bursts is the hardest 
(see figure H), while at the same time being one of the weakest. Since each burst has a 
different duration, the slopes in physical units may not be directly comparable. However, 
the fitted slope in j3 times the duration of the fitted time interval, from column 6 of table [I], 
is a dimensionless parameter (Aj3) that represents the total change in (3, assuming that 
the evolution in f3 is linear (as it is for the majority of the sample). Figure shows that 
the distribution of A/3 has a single, roughly symmetric peak centered on —0.374, with one 
outlier (not shown in the figure). The median absolute deviation width of the distribution 
is ?i>adcv = 0.392, compared with a standard deviation of a = 0.516. This argues that a 
single physical process characterizes the majority of the sample; and again points out that 
hard-to-soft spectral evolution is typical behavior for the high-energy power-law component. 
Physical mechanisms for burst energetics should account for this, possibly via depletion 
of a reservoir of energy that is available for the burst. Otherwise, it may be that when 
the high-energy portion of the emission changes beyond this point, the total emission is 
quenched. 

The linear fit to the power-law indices does not characterize the distribution well 
for many bursts (24 out of 122), indicating that other types of behavior may be present. 
Figure [5| serves as an example of a burst that has strong hard-to-soft spectral evolution 
but where the linear fit is unacceptably poor. The residuals to the fit have considerable 
scatter that is correlated in successive time bins in several places on the figure. It is these 
residual patterns that we are interested in. Two possibilities are easily tested: there may 
be a correlation between the high-energy behavior and intensity within a burst (clearly not 
the case for 3B 911118 in figure ||]), or the high-energy spectrum may be correlated with 
the evolution of the low-energy spectrum. Burst 3B 911118 is an example of this behavior, 
as can be seen in figure ||, where the fitted values of E peilk (representing the low-energy 
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behavior) and (3 have been plotted against each other. 

For the case of correlation between hardness (as measured by the high-energy power-law 
index) and intensity (measured as total count rate in the fitted energy interval: ~ 28 - 
1800 keV), we applied two statistical tests to the data and multiplied their probabilities 
in order to screen for candidates. The tests (described below) are likely to be correlated; 
however, each measures the hardness-intensity correlation differently, so that their product 
combines the best of each. We set the threshold for significance at 10~ 6 for the product, 
so to avoid false positives as much as possible. In both tests, we removed the first-order 
trend in the data by dividing by the linear fit to the power-law indices (which is described 
above). We do this, despite the fact that many bursts don't show a linear trend in the 
high-energy power-law index, since there are a considerable number of bursts that do have 
a significant correlation between (3 and time, while the burst intensity manifestly does 
not: a typical burst will have overlapping regions of both positive (rising portions) and 
negative (falling) correlation with time, so that the whole ensemble of (3 values has no 
correlation. The overall linear trend may be larger than the amplitude of the residuals of 
the fitted linear model (this is the case in figure H), in which case there is no significant 
hardness - intensity correlation as determined by (3 alone. After detrending, the residuals 
may or may not be correlated with intensity. The Spearman rank-order test is relatively 
unequivocal: that is, if the resulting probability is low enough, then the desired correlation 
definitely exists. However, the converse is not true: the test can fail badly since it ignores 
the one-sigma errors in the fitted power-law indices. For this reason, we also have calculated 
the linear correlation coefficient between the detrended values of (3 and intensity, where 
the inverses of the variances on the detrended power-law indices are used to weight their 
contribution ( [Press et al. 199 1 ). For this case, individual, poorly-determined indices that 
are only a few sigma away from being consistent with correlation contribute the same as 
well-determined ones closer to the center of the distribution. In practice, while this kind of 
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test is a poor indicator of whether an observed correlation is statistically significant, it is 
a rough indication of the strength of a correlation under the assumption that a correlation 
definitely exists, so the two statistical tests we've chosen complement each other, to a 
certain extent. Their product selects those bursts that have low probabilities (indicating 
strong correlation) from both tests (assuming that by detrending no significant correlation 
was introduced that was not present in the original data). We have indicated the combined 
probabilities from both tests in the seventh column of table ^ and also indicate the sign 
of the linear correlation coefficient. Since the power-law indices were detrended, a positive 
sign indicates a negative actual correlation; that is, the high-energy behavior is opposite 
that of the burst time history. An example of positive correlation in the detrended values 
of (3 for 4B 960924 is shown in figure || A small number of bursts (9), have significances 
less than 10~ 6 . Of these, 6 are examples of positive correlation. A larger number (24) are 
significant at the 10 -4 level. 

Another possible type of behavior in (3 that is testable with our data is a correlation 
with the low-energy spectral evolution. The most obvious such behavior is the hard-to-soft 
spectral evolution of E peak , discussed by [Ford et al. (19951 - -^pcak is a good measure for 



overall spectral evolution since it marks the peak in the power output of the spectrum per 
log. decade. Of course, -E pea k is not defined for those portions of a burst where (3 > —2; in 
that case, we substitute the break energy of the spectrum instead. In addition, we wish to 
check for higher moments of correlation than is possible with a linear trend of (3 in time, 
which was discussed above, such as the evolution of (3 within individual peaks of a burst. In 
table column 8, we calculate the Spearman rank-order probability that the distribution 
of j3 for a given burst is correlated with the distribution of -Epcak; which stands in here for 
the low-energy behavior. The best example of correlation with E pea ^ is shown in figure || 
which is a plot of the two fitted parameters against one another for 3B 911118. Out of 122 
bursts, 15 bursts have probabilities less than 10~ 3 , indicating correlation, and out of these 
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only 5 have significant hard-to-soft spectral evolution, as measured by how many sigma 
the slope in f3 in table |2|, column 5, deviates from 0. The important point is that, whereas 
hard-to-soft behavior can be demonstrated for large numbers of bursts in the evolution both 
i?pe a k and P, this behavior is generally not correlated between the two. Indeed, hard-to-soft 
evolution of -E pea k within individual peaks of a burst is not typically observed with /3, 
otherwise, far more instances of correlation between the two would have been observed. 



4. Discussion 

In this series of BATSE spectral analysis papers, we have demonstrated several times 
the universal suitability of the 'GRB' spectral form for fitting burst spectra, whether it is 
applied to the total spectrum averaged over the burst ( [Band et al. 1993|) , to time- resolved 



spectroscopy of bright bursts in the SD data in [Ford et al. (1995)] , to joint fits of 



time-averaged spectra of bright bursts with the low-energy discriminator data (preece et 



al. 1996| ; although we see the model break down with low-energy excesses observed in 



15% of GRBs) and now to time-resolved spectroscopy of bursts observed mostly with the 
BATSE LADs. In figure |], we now see that there is evidence of an average high-energy 
power law index that is ~ —2 in a large number of GRBs. In addition, the variance of this 
index over the sample is similar to that obtained by Pendleton ( 1994a]] , using BATSE LAD 
discriminator data. 

Table presents evidence that (3 is not constant for 42 out of 122 bursts in our 
sample. The typical change in {3 over an entire burst, ~ 0.4 (figure 0), is small compared 
with the average value of /3 ~ —2.1. We should consider which of the many emission 
models proposed for GRBs are consistent with these observations. A —2 power law slope 
is evidence for single-particle cooling, from either synchrotron losses or Compton scattering 
( |Blumenthal fc Gould 19TD| ). Typically, one would integrate the energy loss rate over the 
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particle distribution; however, particles that are relatively cool with respect to their large, 
possibly relativistic bulk motion can be treated as monoenergetic in interactions with static 
external particles or fields. Bremsstrahlung losses are another matter. Such scenarios have 
been proposed for bursts of cosmological origin for external shocks, ( [Rees fc Meszaros 1992 



& [Meszaros fc Rees 1993|) as well as for synchrotron shocks QKatz 1994] & |Tavani 1996|) . It 



should be noted that the cooling timescale for most expected processes, especially those 
like synchrotron that involve magnetic fields, is far shorter than observed burst lifetimes 
by many orders of magnitude. In fact, this is a common problem with GRB models: 
an unspecified energy storage mechanisms usually must be invoked in order to extend 
the emission. Relativistic bulk motion, which is necessary to ensure that bursts do not 
degenerate into a pair fireball, can multiply the lab-frame lifetime by the Lorentz factor, 
usually considered to be on the order of 1000. This is not nearly long enough for processes 
such as synchrotron emission whose characteristic timescale may be on the order of 10~ 17 s. 
Clearly, in bursts, there is a reservoir of energy, possibly the protons that carry the bulk of 
the kinetic energy in the blast wave. 

It appears that hard-to-soft spectral evolution predominates over soft-to-hard, as 
observed already in N orris et al. (1986) , Ford et al. (1995) and Band (1997) . In our 



study, the high-energy behavior follows this trend at the greater than the 3a level in 50 
out of 122 cases, while the opposite is true for only 5 bursts at the same significance. This 
is independent of the low-energy behavior; indeed, we have a significant correlation with 
the low-energy behavior in only 15 cases and out of these, 5 have significant hard-to-soft 
spectral evolution. Taken together, we have evidence that the high-energy behavior is very 
much independent of the rest of the spectral evolution of a burst; in 35% of the cases, there 
is hard-to-soft spectral evolution, and no evolution in most of the rest, only 10% of all 
bursts failing the linear fit x 2 test- 
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As seen in figure [|, there is a small group of 'super-soft' bursts characterized by 
Awe <~ —3.0. Along with 4B 970111, which was an extremely bright burst with no 
apparent high-energy power law component (it was fitted with the COMP model), we have 
three such events. Several of these have no detectable emission above ~ 600 keV. This 



behavior is similar to the 'no high-energy' bursts of |Pendleton et al. ( 199771 which were 
shown to be homogeneous in space. Since most of the homogeneous bursts were relatively 
weak, compared with the entire sample, here we must be observing the brightest few of 
that set, rather than 20%, as reported in Pendleton et al. (1997) . There may actually 
be a continuum of burst properties, with these bursts representing the furthest extreme. 
Bursts in this extreme (as well as some portions of other bursts that have very steep 
high-energy power laws) may be an indicator that some emission-limiting phenomenon 
such as pair-plasma attenuation may be at work. Indeed, in many cases, spectra in these 
bursts can be fitted by a spectral form that does not require a high-energy power-law (such 
as the COMP model). This also fits in with the observation that such events are typically 
weaker than average. In the context of shock models of GRBs, several parameters of the 
particle energy distribution determine the resulting spectrum. These may be factors such 
as the shape of the distribution, whether it is a power law, the maximum energy or the 
bulk Lorentz factor. It is very likely that the maximum energy of the accelerated particle 
distribution resulting from the shock could be drawn from an enormous range (out to 
several GeV, at least), depending on the conditions at the shock. Thus, the super-soft 
bursts may be representative of particle distributions that arise from weak shocks, affecting 
the shape or maximum energy in such a way to limit the high-energy emission. 
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5. Summary 

In this study, we have looked in detail at the temporal behavior of the high-energy 
power-law portion of GRB spectra from a sample of 126 bursts selected by either high 
flux or fluence. The average over all fitted spectra of all bursts in the sample for the high 
energy power law index (/?) is ~ —2.12, although fitting a constant, average index to the 
time history of (3 in each burst resulted in unacceptable x 2 values for 34% of the bursts. In 
addition, of those bursts in which f3 is not constant, a large number (100) show hard-to-soft 
spectral evolution, compared with those that have an overall, significant soft-to-hard 
trend. The total change of (3 over the time interval chosen for fitting has a single-peaked 
distribution, centered on —0.37, indicating that theoretical modeling will have to explain 
why most bursts favor this value. In several bursts, the hard-to-soft spectral evolution is 
correlated with similar behavior at lower energies. We also find that some bursts have a 
significant correlation between f3 and the burst time history, or equivalently, instantaneous 
flux. Some bursts in the sample were too soft to be characterized by a high-energy tail, 
while there are intervals in many bursts that have similar behavior, as has been reported by 
Pendleton et al. (1997) . Taken together, these results show that the high-energy spectral 



component has a rich life, independent to a large extent of the behavior of much the rest of 
the spectra. 
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BATSE webserver: |http: / / www.batse.msfc. nasa.gov/| . BATSE work at UCSD is supported 
under NASA contract NAS 8-36081. 
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A. BATSE Large Area Detectors 

The BATSE LADs are a set of eight identical Nal detectors, which are mounted on 
the corners of the CGRO and oriented to ensure maximum all-sky exposure. Perhaps 
the most important feature of the BATSE instrument is its ability to localize a transient 
cosmic source by the comparison of counting rates in the four detectors that directly see it 
( [Pendleton, Briggs fc Meegan 19961) . This is an invaluable aid to spectroscopy, since the 



detector response is a strong function of the source-to-detector axis angle, with differing 
responses at different energies ( jPendleton et al. 1995| ). Thus, without location information, 
the detector response cannot be fully modeled, and spectral model fitting cannot be done 
accurately 

Spectral data from the LADs are compressed to either 128-channel high energy 
resolution background and burst data (HER and HERB datatypes, respectively) or 
16-channel continuous background or medium energy resolution burst data (CONT and 
MER datatypes). The HER background data are typically accumulated over 300 s, while 
the CONT data are always accumulated every 2.048 s. The HERB burst trigger data 
are accumulated in a time-to-spill mode: one spectrum is generated in the time it takes 
to record 64 k counts (in units of 64 ms), currently with a fraction of the last available 
background rate subtracted, to ensure that longer accumulations are taken over periods 
when the burst has returned to background levels. This fraction was zero for roughly the 
first half of the mission, so bright, highly- variable bursts commonly ran out of available 
memory. For the four detectors recording the highest count rates at the time of the trigger, 
there are 128 spectral accumulations, each 128 ms in duration or greater. The lowest seven 
channels of the 128 are at or below the analog lower- level discriminator (LLD), and are 
unusable; the highest few channels suffer from saturation in the pulse amplifier and thus are 
also thrown out. The remaining channels are spaced quasi- logarithmically in energy, falling 
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between approximately 28 keV and 2 MeV, with the exact energy coverage of each channel 
in each detector determined by a channel-to-energy conversion algorithm. It is important 
to note that these energy ranges are quite stable through the mission, due to automatic 
gain control of the PMT voltages. The energy resolution of the LADs was measured on the 
ground to be ~ 20% at 511 keV ( [Horack 1991] ) , and has been quite stable in orbit. 

B. Energy Calibration Methodology 

In order for spectroscopy to be possible with the LAD HERB data, we have had to 
apply a correction to the channel-to-energy conversion algorithm that was developed before 
the launch of the spacecraft. Measurements of several calibration sources at known energies 
resulted in an empirical relationship between channel number and channel energy threshold 
( [Lestrade 199 1|) . The function fitted was essentially linear, with a small non- linear term 
(significant only at low energies), proportional to the square root of the channel number; 
thus there are three fitted parameters. After several bright bursts were observed in orbit, 
it became clear that each detector had a systematic pattern of residuals, localized to the 
low end of the count spectrum. With the assumption (tested below) that these features are 
intrinsic to each detector, and not a function of detector-to-source angle or source intensity, 
we developed a method of calibration using in-orbit data. 

In order to properly calibrate the detectors, we must choose bright objects with 
well-known spectral properties, seen by each detector. Solar flares are generally not usable, 
since they are rarely seen by half the detectors, due to the pointing constraints of the 
spacecrafts solar panels, and their spectra are typically too soft. Earth occultation data 
from the Crab nebula was used by [Pendleton [1994b] to calibrate the 16 channel CONT 
spectra. However, this was not feasible for HER spectra because of telemetry constraints. 
We are left with bursts themselves. Averaged over their entire time history, at least some 
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bursts can be expected to have a fairly smooth spectrum ( [Band et al. 1993 ). Spectral 



features, such as lines, will tend to average out over time and in the LAD data will not 
contribute much overall, due to the moderate energy resolution of the detectors. For 
bright bursts, we can precisely determine the average spectrum from the well-calibrated 
SD spectral data (|Band et al. 1992] ) to use as a constraint in a joint fit with the LAD 



spectral data. The single time-averaged spectrum from the calibration burst is no longer 
available for spectroscopy; however, individual spectra from the burst are still usable for 
our analyses, for two reasons. First, continuum spectral fits are robust, in that they sample 
broad features in the spectrum, rather than the behavior of individual channels. Second, 
the calibration affects only the lowest channels of the spectra, and therefore does not affect 
spectral fitting of the high-energy power-law index, as long as it is determined by counts 
above ~ 150 keV. In the present paper, we needed to obtain a global fit to each spectrum, 
so it was important to calibrate the lower channels as well as possible. 

The general process is iterative: we jointly fit the LAD and SD spectral data for an 
entire outburst interval in a bright burst, using the standard calculation for the LAD data 
energy thresholds. The residuals of the fit to the LAD data are used to determine by how 
much to adjust the energy of each data channel edge in order to bring the count rate closer 
to the model rate. With this new set of edges, a new detector response matrix (DRM) is 
generated to account for the shift in the position of the photopeak with the change in output 
edges and the accompanying change in total response. The photon model is recalculated 
with the new DRM and count rate residuals are again determined. Since the pre-flight 
calibration produces acceptable agreement above ~ 150 keV, we limited the re- calibration 
to energies below 150 keV. We also enforced a fixed lower energy for HERB data channel 
7 (~ 25 keV), to limit the corrections to apply only to channels above the energy of the 
LLD, as this is currently not modeled in the DRM. The freedom of lower-energy edges 
to wander is highly constrained in the joint fit with the SD data, which overlap the LAD 
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energy range and can extend the continuum fit to lower energies by up to 10 keV. Each 
of the edges within the two limits are recalculated in each new cycle until the value of x 2 
for the fit stops decreasing. For each of the eight detectors, one calibration burst yields a 
set of offsets of new edges relative to the original edges, which can then be applied to all 
bursts observed by that detector throughout the mission. We have extensively tested the 
hypothesis that the non-linearities are intrinsic to the detector by examining the residuals 
to spectral fits of several very bright bursts in each detector with the new calibration. We 
have found excellent agreement of the calibration results between bursts, regardless of the 
angle, intensity or hardness of any given event. 
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Fig. 1. — A comparison between the fitted values for (3 from two different spectral models 
(GRB & BPL) for 4B 950403. 

Fig. 2. — A 2-D contour plot showing correlation between -Z? pea k and (3 in the broken power- 
law model. A solid dot indicates the best-fit values for the two parameters. As discussed 
in the text, note that the la contour on this figure is appropriate for only one parameter 
of interest, that is, it represents A% 2 = 1, rather than 2.3, which would contain 68% of the 
joint probability. The data are from the interval 7.680 - 7.808 s of 4B 950403. 

Fig. 3. — Histograms of the log. of the x 2 probability that (3 is a constant (dotted line) or 
exhibits a linear trend (solid line) throughout each burst. The bursts that have probabilities 
consistent with zero (indicating rejection of the model) are not shown for either of the two 
distributions (34 and 20, respectively). 

Fig. 4. — Histogram of the weighted average of (3 for the burst sample (solid line). 
Overplotted is the subset of bursts for which a constant value for (3 resulted in acceptable 
values for x 2 , at the 10~ 4 level of significance (dotted line). 

Fig. 5. — Example of hard-to-soft spectral evolution in (3 for 3B 911118. Both the fitted 
value of (3 with their errors (solid crosses) and the burst count rate history (dotted lines) are 
plotted as a function of time. 

Fig. 6. — Histograms for the distributions of coefficients of correlation between (3 and time 
(solid line) and E pca ^ (dotted line) for the burst sample. 

Fig. 7. — Histogram of the total change in the high-energy power-law index (fitted value of 
dp/dt times the total time interval) for each burst in the sample. To improve the clarity of 
the figure, one outlier is not shown. 

Fig. 8. — The fitted values of E peak plotted against (3 for 3B 911118, illustrating correlation 
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between the time evolution of both hard and soft spectral components in a burst. 

Fig. 9. — Example of positive correlation between fitted values of f3 divided by their linear 
trend ('detrended' f3 - crosses) and the time history of the count rate for 4B 9960924 (the 
count rate has been divided by the energy range 28 — 1800 keV - solid histogram). 



